Bose-Einstein condensation in quantum glasses 
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The role of geometrical frustration in strongly interacting bosonic systems is studied with a com- 
bined numerical and analytical approach. We demonstrate the existence of a novel quantum phase 
featuring both Bose-Einstein condensation and spin-glass behaviour. The differences between such 
a phase and the otherwise insulating "Bose glasses " are elucidated. 



Introduction. Quantum particles moving in a disor- 
dered environment exhibit a plethora of non-trivial phe- 
nomena. The competition between disorder and quan- 
tum fluctuations has been the subject of vast literature 
[l|, in past years, with a renewed interest following 
from the exciting frontiers opened by the experimental 
research with cold-atoms 0, Q . One of the most striking 
features resulting from the presence of a disordered ex- 
ternal potential is the appearance of localized states 1 1 . 
Localization happens both for fermions and bosons 
but in the latter case one has to introduce repulsive inter- 
actions to prevent condensation of particles in the lowest 
energy state. This results in the existence of an insu- 
lating phase called "Bose glass", characterized by a finite 
compressibility and gapless density excitations in sharp 
contrast to the Mott insulating phase [3, [1] . 

On the other hand, latest research stimulated by the 
discovery of a supersolid phase of Helium has brought 
to the theoretical foresight of a "superglass" phase [l,13|i 
corroborated by recent experimental evidence Q , where 
a metastable amorphous solid features both condensa- 
tion and superfluidity, in absence of any random exter- 
nal potential. The apparent irreconcilability, between the 
current picture of insulating "Bose glasses" and the emer- 
gence of this novel phase of matter, calls for a moment 
of thought. Although it has been recently demonstrated 
that attractively interacting lattice bosons can overcome 
the localization induced by an external random potential 
and feature a coexistence of superfluidity and amorphous 
order , a general understanding of the physics of Bose- 
Einstein condensation in quantum glasses and in presence 
of purely repulsive interactions is still in order. In partic- 
ular, we wonder what could be the possible microscopic 
mechanism leading to super-glassines and if the external 
disorder, current paradigm in the description of quantum 
glasses, could be replaced by some other mechanism. 

In this Letter we show that geometrical frustration 
is the missing ingredient. Geometrical frustration is a 
well recognised feature of disordered phases in which the 
translational symmetry is not explicitly broken by any 
external potential. Examples are spin liquids phases of 



frustrated magnets [lO|, valence-bond glasses and 
the order-by-disorder mechanism inducing supersolidity 
on frustrated lattices Another prominent manifes- 

tation of frustration is the presence of a large number 
of metastable states that constitutes the fingerprint of 
spin-glasses. When quantum fluctuations and geomet- 
rical frustration meet, their interplay raises nontrivial 
questions on the possible realisation of relevant phases 
of matter. Most pertinently to our purposes: can quan- 
tum fluctuations stabilise a superglass phase in a self- 
disordered environment induced by geometrical frustra- 
tion? Hereby we answer this question demonstrating 
that repulsively interacting bosons can feature a low- 
temperature phase characterised both by spin-glass or- 
der and Bose-Einstein condensation. Such a frustration 
induced superglass sheds light onto a novel mechanism 
for glass formation in bosonic systems noticeably differ- 
ent from the localization effects leading to "Bose glass" 
insulators and paving the way to a better understanding 
of this new phase of the matter. 

Model. Strongly interacting bosons on a lattice can 
be conveniently described by means of the extended Hub- 
bard Hamiltonian, namely 
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where bl{bi) creates (destroys) a hard-core boson on site 
i, Hi = b\bi is the site-density and the summations over 
the indexes {i,j) are extended to nearest-neighbouring 
vertices of a given lattice with L sites. In the following 
we will set t = 1, i.e. we will measure all energies in units 
of t. In this work, to capture the essential physics of the 
problem in exam, we adopt a minimal and transparent 
strategy to induce geometrical frustration in the solid. 
We therefore consider the set of all possible graphs of L 
sites, such that each site is connected to exactly z = 3 
other sites, and give the same probability to each graph 
in this set. We will discuss average properties over this 
ensemble of random graphs in the thermodynamic limit 
L oo. The motivations for this choice are the follow- 
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Figure 1: Edwards- Anderson order parameter (top) and con- 
densate fraction pc/p (middle) as functions of V at half- 
filling, computed via the cavity method at different values 
of p. In the middle panel Pc/p as obtained by SGF at 
/3 = 5 is reported. (Bottom) Scaled spin-glass susceptibility 
XsG ~ Xsg/L"'^^ reported as a function of V; standard finite- 
size-scaling arguments [S^l show that the different curves must 
intersect at the spin-glass transition. 



ing: i) On a square lattice, model ([T]) is known to produce 
a solid insulating phase at high enough density, where 
the particles are arranged in a checkerboard pattern [31 ■ 
This is due to the fact that all loops have even length. 
On the contrary, typical random graphs are character- 
ized by loops of even or odd length; in the classical case 
t = 0, this frustrates the solid phase enough to produce a 
thermodynamically stable glass phase at high density 14 1 . 
a) Typical random graphs have the important property 
that they are locally isomorphic to trees, since the size 
of the loops scales as InL for large L: indeed, this is a 
consistent way of defining Bethe lattices without bound- 
ary This locally tree-like structure allows to solve 
the model exactly, at least in the liquid phase, by means 
of the cavity method 13, 3|. Hi) These lattices are quite 
different from square lattices. Yet, it has been shown in 



for the model defined on a random graph and on a square 
lattice 15, [11]. Hence, we believe that it is possible to 
find a model similar to Eq. ([1]), defined on a square lattice 
but with slightly more complicated interactions (proba- 
bly involving many-body terms) that will show the same 
qualitative behaviour of the model investigated here. 

Methods. The stochastic sampling of the quantum 
partition function Z = Tre^^'^ at finite temperature 
T = 1/(3 can be conveniently exploited to obtain nu- 
merically exact properties of a generic bosonic Hamilto- 
nian such as ([T|). Quantum Monte Carlo (QMC) schemes 
based on the original Worm algorithm idea [l9| have been 
recently extended to Canonical ensemble simulations [H, 
[2l| . These methods offer an efficient scheme based on the 
sampling of the configuration space spanned by the ex- 
tended partition function Z„(t) = Tre^f'^^^^^M^e^^'", 
where is a suitable worm operator determining an 
imaginary time discontinuity in the world-lines. We have 
chosen the worm operator introduced in 2l|, which is a 
linear superposition of n-body Green functions, avoiding 
the complications arising in [20[ where the commutability 
of the worm operator with the non-diagonal part of the 
Hamiltonian is required. Full details of the Stochastic 



the classical case, and for some more complicated interac- 
tions, that the phase diagram is qualitatively very similar 



Green Function (SGF) method arc described in Ref. [2l|, 
we only stress here that access to exact equal-time ther- 
mal averages of ri-body Green functions is granted as 
well as to thermal averages of imaginary time correla- 
tion functions of local, i.e. diagonal in the occupation 
numbers representation, quantum operators. 

A different and complementary approach to models de- 
fined on random lattices consists in solving them exactly 
in the thermodynamic limit L — > oo, by means of the cav- 
ity method [l3| ■ Since local observables are self-averaging 
in this limit, this results in automatically taking into ac- 
count the average over the different realisations of the 
random graphs. For bosonic systems, the cavity method 
allows to reduce the solution of the model to the prob- 
lem of finding the fixed point of a functional equation for 
the local effective action, in a similar spirit to bosonic 
DMFT. All the details of the computation have been dis- 
cussed in [3l , where it has been shown that the method 
allows to compute the average of all the relevant observ- 
ables. However, in the simplest version discussed in [ll |. 
the cavity method can only describe homogeneous pure 
phases such as the low-density liquid. In order to de- 
scribe exactly the high density glassy phase, where many 
different inhomogeneous states coexist, one has to intro- 
duce a generalization of the simplest cavity method which 
goes under the name of replica symmetry breaking (RSB). 
Unfortunately, this is already a difficult task for classical 
models, in particular in spin-glass like phases [l^ . Hence, 
in this paper we describe the glassy phase using the sim- 
plest version of the method, the so-called replica sym- 
metric (RS) one. This yields an approximate description 
of the glassy phase which we expect to be qualitatively 
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Figure 2: Finite temperature phase diagram at half-filling. 

correct. To summarize, in the low-density liquid phase 
we can compute averages numerically with SGF and ana- 
lytically with the cavity method, and we obtain a perfect 
agreement between the two results. In the glassy phase, 
the RS cavity method is only approximate, an exact so- 
lution for L — > oo requiring the introduction of RSB. 
On the other hand, SGF is limited for large L by the 
unavoidable divergence of equilibration times due to the 
glassy nature of the system. Still, we find a good agree- 
ment between the result of SGF for fairly large L, where 
the system can still be equilibrated, and the RS cavity 
method for L — > oo, making us confident that the qual- 
itative and quantitative picture of the glassy phase we 
obtained here is fully consistent. Moreover, we solved 
the model at the simplest (one-step) RSB level in some 
selected state points and we found a very small quanti- 
tative difference with the RS solution. 

Results. The presence of off-diagonal long range or- 
der can be conveniently detected by considering the large 
separation limit of the one-body density matrix, i.e. the 
condensate density reads 
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where the square brackets indicate a quantum and ther- 
mal average and the bar indicates averages over inde- 
pendent realizations of the random graphs. The cav- 
ity method works in the grand-canonical ensemble and 
gives direct access to the average of 5, while canonical 
ensemble simulations done with SGF give easy access to 
the one-body density matrix. On the other hand, spin- 
glass order is signaled by the breaking of translational 
invariance, namely {tii) ^ X]i=i("-i) ~ P- Introduc- 
ing 5ni = {jLi — p), the on-site deviation from the av- 
erage density, spin-glass order can be quantified by the 
Edwards- Anderson order parameter 
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which can be easily computed by the cavity method, or 
by the divergence of the spin-glass susceptibility 
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which is more easily accessible in SGF. It is possible to 
show [221 that Xsg is the susceptibility naturally associ- 
ated to the order parameter (Jea, because it can be de- 
fined as the derivative of ^ea with respect to an external 
field coupled to the order parameter itself (as in standard 
critical phenomena). 

At half-filling factor p = 1/2, the condensate fraction, 
the Ed wards- Anderson order parameter, and the scaled 
spin-glass susceptibility are shown in Fig. [TJ In middle 
panel we compare the values of condensate fraction ob- 
tained via the cavity method and via SGF in a linear 
extrapolation to L ^ oo. The very good coincidence 
of these results supports our conjecture that the approxi- 
mate RS description of the glass phase we adopted here is 
quantitatively and qualitatively accurate. At the lowest 
temperature, the system becomes a glass around V ^ 2.7 
while it still displays BEG; the condensate fraction only 
vanishes at V ~ 3.5 inside the glass phase. This clearly 
establishes the existence of a zero-temperature supcrglass 
phase in the region 2.7 F ^ 3.5. Note additionally that 
both transitions are of second order, hence the conden- 
sate fraction is a continuous function; since the latter 
stays finite on approaching the spin-glass transition from 
the liquid side (where the cavity method gives the exact 
solution), it must also be finite on the glass side just after 
the transition. In Fig.[2]we report the finite temperature 
phase diagram of the model at half-filling. It is defined by 
two lines: the first separates the non-condensed ((6)=0) 
from the BEG ((&) ^ 0) phase, the second separates the 
glassy ((7ea 0) from the liquid ((7ea = 0) phase. The 
intersection between these two lines determines the exis- 
tence of four different phases (normal liquid, superfluid, 
normal glass, superglass). 

Ground-state degeneracy. Geometrical frustration in- 
duces the existence of a highly degenerate set of ground- 
states, each of them characterized by a different average 
on-site density, which is absent in glassy phases induced 
by localization in disordered external potentials such as 
the Bose glass. To demonstrate this peculiar feature, it is 
instructive to consider a variational wave-function explic- 
itly breaking the translational symmetry of the lattice 
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where the variational parameters are explicitly site- 
dependent and tend to (dis)-favour the occupation of a 
given site. In the spin-glass phase of the bosons, the opti- 
mal set of the variational parameters is highly dependent 
on the initial conditions associated with the a, , whereas 
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Figure 3: Variational expectation values of the site density 
for different sets of the optimized parameters at half-filling 
density for L = 80 and V = 4. 



all the variational states, even with different parameters, 
have almost degenerate variational energies. Each set of 
optimized variational parameters is then representative 
of one of the many degenerate ground-states of H. As 
an example, we show in Fig. [3] the variational expecta- 
tion values of the site densities for two different solu- 
tions resulting from the minimization of the variational 
energy with the SRH method [1^, a robust stochastic 
variant of the Newton Method. We further checked, us- 
ing the zero-temperature Green Function Monte Carlo 
method that if one applies the imaginary-time evo- 
lution \(j)a) = exp(— r'H)|vE'Q) to one of these states, the 
density profile remains amorphous for a time r that is 
divergent with the size of the system. 

Conclusions. The aim of this Letter was to show the 
existence of a stable superglass phase in a lattice model of 
geometrically frustrated bosons, in absence of quenched 
disorder in the Hamiltonian. This has been done by 
combining the analytical solution of the model via the 
quantum cavity method and numerical simulations via 
Quantum Monte Carlo. The glass phase we found is 
very different from the usual Bose glass, since the latter 
is driven by localization effects in presence of an exter- 
nal disorder and is then insulating, while the former is 
driven by self-induced frustration on a disordered lattice 
and displays Bose-Einstein condensation. This results 
in a coexistence of a large number of degenerate amor- 
phous ground states, whose existence we showed by a 
variational argument corroborated by QMC. We expect, 
by analogy with the classical case [15|, that the glassy 
phase found here will exist also on regular finite dimen- 
sional lattices, provided the interactions are modified to 
induce sufficient geometrical frustration. In that case, its 
properties should be very similar to the one showed by 
metastable superglasses observed both in numerical sim- 
ulations 0] and experiments Q on Helium 4. The main 
difference is that, due to the randomness of the under- 
lying lattice, the superglass studied here is a truly stable 
equilibrium state, allowing for a much more precise char- 



acterization of its properties. 
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Note added. After this work was completed, we be- 
came aware of the paper (2^ where related results have 
been obtained. 
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